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Abstract 

A  procedure  is  developed  to  analyze  uncertain 
systems  having  both  time-varying  (TV)  and  lin¬ 
ear,  time-invariant  (LTI)  uncertainty.  The  prob¬ 
lem  is  formulated  as  a  two  step  procedure  involv¬ 
ing  the  solution  of  two  sets  of  LMIs  using  both 
constant  and  frequency  dependent  scales  to  ac¬ 
count  for  the  different  types  of  uncertainty.  In 
the  first  step,  the  scales  corresponding  to  the 
TV  parameters  are  constrained  to  be  constant 
over  a  set  of  frequencies  and  both  constant  and 
frequency-dependent  scales  are  computed  to  min¬ 
imize  the  peak  value  of  the  robustness  bound  for 
the  chosen  frequencies.  In  the  second  step,  the 
constant  scales  are  held  fixed  and  the  frequency- 
dependent  scales  corresponding  to  the  LTI  un¬ 
certainty  are  computed  to  minimize  the  robust¬ 
ness  bound  at  each  frequency.  Numerical  ex¬ 
amples  are  given  to  demonstrate  the  procedure, 
which  gives  a  less  conservative  result  than  pre¬ 
viously  shown  and  thus  allows  for  more  accu¬ 
rate  analysis  of  systems  with  TV  and  LTI  un¬ 
certainty. 

1  Introduction 

Control  law  design  for  systems  operating  over  a 
wide  range  of  conditions  is  traditionally  a  heuris¬ 
tic  process  consisting  of  several  steps.  First,  con¬ 
trol  laws  are  synthesized  for  a  representative  set 
of  operating  conditions  by  fixing  the  system’s  pa¬ 


rameters  in  each  design.  The  resulting  set  of  con¬ 
trol  laws  is  then  implemented  in  the  parameter- 
varying  system  by  scheduling  with  operating  con¬ 
dition.  This  gain  scheduling  can  consist  of  either 
switching  between  fixed  control  laws  or  transi¬ 
tioning  gradually  from  one  to  another  using  in¬ 
terpolation.  Stability  of  the  closed-loop  system 
is  achieved  by  sufficient  robustness  margins  at 
each  design  point  and  can  be  evaluated  through 
exhaustive  simulation,  although  there  is  no  guar¬ 
antee  of  stability  for  all  operating  conditions  and 
time-consuming  redesign  is  often  necessary.  Fur¬ 
thermore,  as  the  complexity  of  the  system  in¬ 
creases  or  as  more  sophisticated  control  law  de¬ 
sign  techniques  are  used,  these  steps  can  become 
difficult,  as  the  resulting  control  laws  can  become 
more  difficult  to  gain  schedule  because  of  their 
increased  or  even  varying  dynamic  order. 

To  overcome  these  obstacles,  Packard  [7]  de¬ 
veloped  a  technique  to  address  the  problem  of 
gain  scheduling  controllers  by  representing  the 
plant  as  a  linear  parameter-varying  (LPV)  sys¬ 
tem.  Here,  the  parameter-varying  system  is  rep¬ 
resented  as  the  linear  fractional  transformation 
of  a  nominal  system  with  the  time-varying  pa¬ 
rameters.  The  controller  is  then  parameterized 
using  the  same  set  of  time-varying  parameters 
so  that  the  resulting  controller  is  scheduled  a 
priori  with  operating  condition  and  hence  the 
problem  of  gain  scheduling  a  set  of  LTI  con¬ 
trollers  can  be  eliminated.  The  TV  parameters 
in  the  plant  and  the  controller  parameterization 
are  combined  and  the  problem  is  reformulated 
into  the  paradigm  of  a  nominal  plant  with  linear 
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fractional  uncertainty.  This  formulation  allows 
small  gain  synthesis  techniques  such  as  1-Lcc,  or 
synthesis  to  be  used.  However,  when  this  LPV 
structure  is  combined  with  the  traditional  LTI 
uncertainties  such  as  actuator  uncertainty,  sen¬ 
sor  uncertainty,  unmodeled  dynamics,  or  a  fic¬ 
titious  performance  block  the  resulting  system 
has  both  TV  and  LTI  uncertainty  and  the  tradi¬ 
tional  'Hoo  or  ^-synthesis  framework  cannot  ac¬ 
count  for  both  simultaneously.  Hence,  there  is 
a  need  for  a  nonconservative  method  of  analyz¬ 
ing  such  systems.  While  this  framework  was  the 
main  motivation  behind  developing  a  less  conser¬ 
vative  analysis  technique,  there  are  other  types 
of  systems  with  both  TV  and  LTI  uncertainty. 

The  remainder  of  this  paper  is  organized  in 
four  sections.  Necessary  background  theory  is 
discussed  briefly  in  Section  2.  In  Section  3  we 
develop  a  procedure  to  compute  a  robustness 
bound  for  systems  with  both  TV  and  LTI  un¬ 
certainty  that  is  less  conservative  than  currently 
used  methods.  While  not  presented  here,  the 
analysis  procedure  developed  in  Section  3  can  be 
used  with  a  synthesis  procedure  similar  to  DK D 
iteration  presented  in  [9]  to  synthesize  controllers 
that  are  less  conservative  and  provide  better  per¬ 
formance  than  techniques  currently  being  used. 
This  analysis  procedure  is  demonstrated  on  some 
numerical  examples  in  Section  4.  These  exam¬ 
ples  show  that  the  new  method  gives  a  better 
robustness  bound  than  other  techniques.  In  Sec¬ 
tion  5  we  conclude  with  a  brief  summary  and  a 
discussion  of  some  lessons  learned. 

2  Robustness  Analysis 

In  this  section  we  recall  some  basic  results  in  ro¬ 
bustness  theory.  Consider  a  stable  transfer  func¬ 
tion  G(s)  in  a  feedback  interconnection  with  an 
uncertain  matrix  A  G  A,  where  A  has  some 
block  structure  and  o'(A)  <  7"^  for  all  A  G  A. 
The  primary  objective  of  robustness  analysis  is 
to  determine  the  largest  A  for  which  the  sta¬ 


bility  of  the  closed-loop  system  is  guaranteed. 
From  the  small  gain  theorem,  we  know  that  the 
closed-loop  system  is  asymptotically  stable  if 

<  7, 

for  all  w  G  IR.  Thus,  the  maximum  singular 
value,  a(G{jLo)),  gives  a  frequency-dependent  ro¬ 
bustness  bound,  the  peak  of  which  defines  a  lower 
bound  on  the  smallest  destabilizing  A. 

The  above  robustness  measure  can  be  quite 
conservative.  Fortunately,  the  conservatism  of 
this  robustness  bound  can  be  reduced  by  choos¬ 
ing  invertible  scalings  that  commute  with  the 
uncertain  matrix  A.  Let  H  G  X>,  where  V  is 
the  set  of  invertible  scalings  that  commute  with 
every  A  G  A.  Note  that  if  D  G  and  A  G  A, 
then  DAD~^  =  DD~^A  =  A.  In  this  case,  the 
closed-loop  system  is  stable  if  there  exists  D  ^  T> 
such  that 

a{DG{ju})D-^)  <  7 

for  all  w  G  IR. 

Now  define, 

7(0;)  =  a{DG{ju})D-^), 

and 

7  =  sup  7(0;)  =  sup  a-{DG{jw)D~^), 

U)  U) 

so  that  7  is  the  peak  of  7(ct;),  the  scaled  maxi¬ 
mum  singular  value.  Then,  7“^  is  a  lower  bound 
on  the  size  of  the  smallest  destabilizing  A.  Specif¬ 
ically,  stability  is  guaranteed  for  all  A  satisfying 
^(A)  <  7~^,  that  is,  there  is  no  A  satisfying 
a(A)  <  such  that  det(s/  —  (j(_;a;)A)  is  sin¬ 
gular  for  some  w  G  IR.  Then,  the  scalings  can 
be  used  to  minimize  7,  or  equivalently,  to  max¬ 
imize  the  size  of  A  for  which  stability  can  be 
guaranteed. 

Suppose  A  G  A  represents  LTI  uncertainty. 
In  this  case,  the  scaling  matrix  Da  €  D  can  be 
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chosen  to  be  a  function  of  frequency  so  that  the  that  minimize  the  peak  of  the  robustness  bound 
peak  of  the  robustness  bound  is  given  by  given  by 

7  =  sup  d-{DA{i^)G{jLo)D^^{u)).  q  =  sup  a  (^DG{ju>)D~^^ .  (1) 


Hence  the  robustness  analysis  problem  is  to  find 
the  frequency-dependant  scale  D^{u)  that  min¬ 
imizes  the  robustness  bound 

7(0;)  =  a{DAiLo)G{juj)D-^\u:)) 

pointwise  at  each  frequency.  This  can  be  com¬ 
puted  as  the  solution  to  a  convex  optimization 
problem  as  in  [5]. 


Furthermore,  since  the  frequency  dependent  scales 
that  minimize  7  are  only  unique  at  the  frequency 
at  which  7  occurs,  additional  analysis  informa¬ 
tion  can  be  obtained  at  other  frequencies  by  find¬ 
ing  the  frequency-dependent  scales  Da{(^)  that 
minimize 

7H  =  5  {DG(jw)b-')  ,  (2) 


In  contrast,  suppose  0  e  A  represents  TV 
uncertainty.  In  this  case,  the  scaling  Dg  must  be 
fixed  for  all  frequency.  In  this  case  the  peak  of 
the  robustness  bound  is  given  by 

7  =  sup  a(DgG{juj)Dj^). 

Hence  the  analysis  problem  is  to  find  the  con¬ 
stant  scale  Do  that  minimizes  7,  the  peak  of  the 
robustness  bound  7(0;)  =  a{DgG{jw)Dj^)  over 
frequency.  This  problem  can  be  formulated  as  a 
convex  problem  in  the  state  space  as  in  [3].  Al¬ 
ternatively,  if  the  plant  has  the  LPV  form  in  [7] 
the  analysis  can  be  done  as  a  step  of  the  con¬ 
troller  synthesis  using  the  method  presented  in 
[2]. 

In  many  cases,  the  uncertain  matrix  A  G  A 
may  represent  uncertainty  having  both  LTI  and 
TV  blocks.  In  this  case,  the  augmented  uncer¬ 
tainty  A  is 


A  = 


e  0  ■ 

0  A  • 


Invertible  scalings  having  the  complementary  block 
structure  can  be  defined  as 


D  = 


Do  0 
0  Da  {oj) 


The  analysis  problem  is  then  to  find  the  constant 
and  frequency  dependent  scales  Dg  and  Da{i-o) 


at  each  a;  G  IH.  The  method  to  obtain  the  scales 
for  this  case  is  not  as  straightforward  since  some 
of  the  scales  must  be  constant  over  all  frequen¬ 
cies  while  others  can  vary  with  frequency.  The 
procedures  described  for  the  case  of  only  TV  un¬ 
certainty  could  be  applied  here  by  treating  LTI 
uncertainty  as  TV;  however,  this  could  be  con¬ 
servative.  In  order  to  minimize  the  conservatism, 
frequency-dependent  scales  must  be  used  for  the 
LTI  blocks. 

Some  of  the  conservatism  introduced  by  treat¬ 
ing  all  of  the  uncertainty  as  TV  and  using  a 
method  such  as  in  [3]  to  obtain  all  constant  scales 
can  be  reduced  by  performing  a  second  step  which 
scales  the  LTI  blocks  with  frequency-dependent 
scales.  This  3Ds  procedure  is  summarized  be¬ 
low.  Using  a  state  space  formulation,  find  the 
constant  scale  Dg  and  the  constant  scale  Dac 
that  minimize 


sup  (T 

W 


Dg  0 
0  Dac 


GM 


0 


Then,  at  each  frequency  find  the  frequency  de¬ 
pendent  scales,  Z)a(u;),  minimizing, 

7(0;)  =  ^  {dG{jc,)D-^)  . 

Finally,  the  peak  value  of  this  robustness  bound 
gives  7.  This  technique  has  the  advantage  that 
it  can  be  used  to  analyze  any  system  containing 
TV  and  LTI  uncertainties. 


Similarly,  in  [9],  the  method  of  [2]  was  aug¬ 
mented  with  a  second  step  which  scales  the  LTI 
blocks  with  frequency-dependent  scales,  thus  re¬ 
ducing  the  conservatism.  This  DKD  technique 
presented  in  [9]  can  be  summarized  as  follows. 
First,  let  G(5)  =  Fi{P{s),  K{s))  where  P{s)  is 
an  LPV  system  and  K{s)  is  an  LPV  controller, 
both  in  linear  fractional  form  as  described  in 
[7].  Using  the  technique  of  [2],  find  the  constant 
scale,  Dg,  and  the  controller  minimizing 


sup  a 


Dg  0 
0  I 


G{joj) 


scale  must  be  chosen  a  priori.  If  too  low  an  or¬ 
der  is  specified  the  results  could  be  arbitrarily 
conservative.  On  the  other  hand,  the  computa¬ 
tional  time  becomes  prohibitive  for  systems  of 
even  moderate  order.  At  present,  this  technique 
doesn’t  seem  practical  and  was  not  used  for  com¬ 
parison  in  our  numerical  work. 


3  Analysis  with  Varying  and 
Uncertain  Parameters 


Then,  at  each  frequency  find  the  frequency  de¬ 
pendent  scales,  minimizing, 

7(0;)  =  O'  (^DG{ju))D~^'j  . 

Finally,  7  is  given  by  the  peak  value  of  the  ro¬ 
bustness  bound  7(0;).  Note  that  this  procedure 
is  usually  used  in  an  iterative  manner  and  works 
well  for  synthesis  for  the  class  of  problems  de¬ 
scribed  in  [7].  However,  the  analysis  is  a  byprod¬ 
uct  of  synthesis,  which  is  a  rather  limiting  con¬ 
straint  as  often  it  is  desirable  to  analyze  a  con¬ 
troller’s  robustness  to  uncertainties  that  weren’t 
accounted  for  during  synthesis.  Additionally, 
there  are  systems  that  contain  TV  uncertainties 
for  which  the  method  in  [2]  doesn’t  apply. 

Unfortunately,  both  of  the  above  methods 
can  introduce  significant  conservatism,  since  there 
is  no  guarantee  that  the  constant  scale  found 
in  the  first  step  is  the  best  constant  scale  once 
the  system  is  augmented  with  the  frequency- 
dependent  scales  obtained  in  the  second  step. 
This  is  easy  to  verify  by  computing  the  second 
step  using  different  constant  scales.  In  [1]  an  al¬ 
ternative  approach  is  given,  where  the  frequency- 
dependent  scale  is  parameterized  as  a  rational 
function  and  the  optimization  problem  is  refor¬ 
mulated  using  results  such  as  those  in  [4,  8]. 
This  procedure  can  be  very  computationally  in¬ 
tensive  even  for  problems  of  moderate  size,  and 
the  dynamic  order  of  the  frequency-dependent 


In  this  section  we  describe  a  new  procedure  for 
computing  the  optimal  mixed,  constant  and  fre¬ 
quency  dependent  scales  for  systems  with  both 
TV  and  LTI  uncertainty,  that  is,  we  describe  a 
procedure  that  gives  the  least  conservative  re¬ 
sults  possible  for  the  analysis  problem  described 
by  equations  (1)  and  (2).  The  problem  is  for¬ 
mulated  as  a  two  step  procedure  involving  the 
solution  of  two  sets  of  LMIs.  In  the  first  step, 
the  scales  corresponding  to  the  TV  parameters 
are  constrained  to  be  constant  over  a  set  of  fre¬ 
quencies.  Then,  both  constant  and  frequency- 
dependent  scales  that  minimize  the  peak  of  the 
scaled  maximum  singular  value  are  computed 
for  the  chosen  frequencies.  In  the  second  step, 
the  constant  scales  are  held  equal  to  the  value 
found  in  the  first  step  and  the  optimal  frequency- 
dependent  scales  corresponding  to  the  LTI  un¬ 
certainty  are  computed  to  minimize  the  scaled 
maximum  singular  value  at  each  frequency.  It 
should  be  noted  that  the  peak  of  the  robust¬ 
ness  bound  cannot  be  reduced  in  the  second  step. 
Thus,  for  analysis  purposes  where  only  the  peak 
value  of  the  robustness  bound  is  needed,  only 
step  one  needs  to  be  completed. 

Using  the  definition  of  the  maximum  singu¬ 
lar  value,  it  follows  that  there  exists  an  invert¬ 
ible  scaling  D  such  that  a[DG{ju)D~^)  <  7  if 
and  only  if  there  exists  a  positive  definite  scal¬ 
ing  Q  such  that  G*{jlo)QG{joj)  <  where 
Q  E  V  and  D  =  Thus,  the  system  is 
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where 


stable  for  all  A  G  A  satisfying  (t(A)  <  7  Mf 
G*{joj)QG{ju)  <  7^(5,  where  Q  eV. 


Now,  the  analysis  problem  for  systems  with 
both  TV  and  LTI  uncertainties  described  in  equa¬ 
tions  (1)  and  (2)  can  be  restated.  Let  the  scaling 
matrix  Q  be  partitioned  as 


Qe 
0 


where  Q$  is  the  constant  scaling  matrix  corre¬ 
sponding  to  the  TV  parameter  and  is  the 

frequency-dependent  scaling  matrix  that  corre¬ 
sponds  to  the  LTI  uncertainty.  Then,  for  mixed 
TV  and  LTI  uncertainty,  the  analysis  problem  is 
to  find  Qg  and  Q^{u>)  in 


G*{ju:)QG{ju:)  <  ^\lo)Q 


that  minimizes  the  peak  of  7(0;)  for  all  cu  G  IR 
and  minimizes  7(0;)  at  all  other  frequencies.  The 
two  step  procedure  used  to  perform  this  analysis 
can  be  summarized  as  follows.  First  find  Qg  and 
(5a(a^)  minimizing  7  such  that 


G’(]^)QG{]oj)  <  7^(3 


for  all  a;  G  IR.  This  gives  the  optimal  constant 
scale,  Dg  =  Qg^^^.  Then,  at  each  frequency  lo 
find  a  new  Qi\{uj)  that  minimizes  7(a;)  at  that 
frequency  such  that 


G-(,uj)QG{jui)  < 

This  gives  the  optimal  frequency-dependent  scale, 

These  two  steps  are  refered  to  as  the  ^DD’’ 
procedure  and  are  implemented  using  LMIs  as 
follows. 


1.  Choose  a  set  of  frequencies  LOi  and  evalu¬ 
ate  the  system  frequency  response  G(juji) 
at  those  frequencies.  Form  the  system  of 
LMIs 


Qe 

QAi^Oi) 

G*{jUi)Q,Gijcoi) 


>  0, 

>  0,  f  =  1, . . .  ,n, 

^  '1  Qi^  ^  l,...,n. 


Qe  0 

0  (5  a  (wi) 


is  the  scaling  matrix  for  the  fth  frequency. 
Then,  find  the  scalings  Qi,  i  =  l,...,n, 
that  minimize  7  subject  to  the  system  of 
LMIs. 


2.  For  each  frequency,  form  the  system  of  LMIs 

QAi^^i)  >  0, 

G*{ju;i)QiG{jUi)  < 

using  the  value  of  Qg  found  in  the  previous 
step.  Then,  at  each  frequency  find  the 
scaling  QA(<^i)  that  minimizes  7^  subject 
to  the  system  of  LMIs. 


The  first  set  of  LMIs  minimizes  the  peak  of 
the  robustness  bound  over  the  set  of  chosen  fre¬ 
quencies  to  yield  the  suboptimal  value  of  the 
constant  scaling  Qg.  As  the  number  of  frequen¬ 
cies  increases,  Qg  approaches  its  optimal  value. 
The  second  set  of  LMIs  minimizes  the  robustness 
bound  at  each  frequency  u>i  to  yield  the  optimal 
values  of  the  frequency-dependent  scales  QA(<^i)- 


4  Numerical  Examples 

In  this  section  we  show  two  numerical  examples 
using  the  procedure  presented  in  the  previous 
section.  The  first  example  is  used  to  demon¬ 
strate  the  procedure,  while  the  second  exam¬ 
ple  shows  that  the  procedure  is  less  conservative 
than  other  methods  used  to  analyze  this  class  of 
problems. 


4.1  Example  1 
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In  this  example  we  use  a  simple,  low  order  model 
to  demonstrate  the  analysis  procedure  described 


in  Section  3.  The  analysis  was  performed  on  a 
weighted,  closed-loop  LPV  model  of  the  VISTA 
F-16  short  period  dynamics  in  linear  fractional 
form  with  time-varying  parameter  6  and  time- 
invariant  uncertainty  A„„c-  The  LPV  model  is 
valid  for  an  altitude  of  1000  feet  and  a  Mach 
number  between  Mmin  —  0.35  and  M^ax  =  0.65. 
The  controller  minimizes  an  'Hoo  norm  for  the 
linear  model  with  Mach  at  its  nominal  value  of 
0.5.  The  state  space  realization  of  the  model  is 
given  in  the  Appendix. 


The  parameter  0  is  a  2  x  2  time- varying  block 
containing  a  normalized  Mach  number,  SM,  on 
its  diagonal,  so  that 


6M  0 
0  5M  ’ 


where  5M{t)  G  [—1  1]  corresponds  to  a  Mach 
number  between  Mmin  and  Mmax-  As  the  Mach 
number  varies,  the  short  period  model  changes 
to  represent  the  short  period  dynamics  at  dif¬ 
ferent  flight  conditions.  The  uncertainty  Au„c  is 
scalar  and  is  augmented  with  a  full  3x3  ficti¬ 
tious  LTI  performance  block,  Aper/-  Thus,  the 
augmented  uncertainty  A  is 


A  = 


■  e  0 

0  Aunc 

0  0 


0 

0 

^perf 


The  invertible  scalings  associated  with  this  block 
structure  can  be  defined  as 


Db  0  0 

0  0 
0  0  dA2{<^)I 


where  Dg  is  a.  2  x  2  constant  matrix  and  d^i  (a;) 
and  dA2i<^)  are  frequency-dependent  scalars.  By 
letting. 


the  analysis  problem  is  in  the  form  of  equation 
(1).  Then,  the  DD  analysis  technique  presented 
in  Section  3  can  be  directly  applied  to  obtain 
the  suboptimal  scalings  De  and  Da{oj).  In  the 
first  step  of  the  analysis  constant  and  frequency- 
dependent  scales  that  minimize  the  peak  of  the 
robustness  bound  are  found  simultaneously.  This 
produces  the  suboptimal  constant  scale.  Then, 
in  the  second  step,  using  the  suboptimal  con¬ 
stant  scale,  new  frequency-dependent  scales  are 
found  that  minimize  the  robustness  bound  at 
each  frequency.  This  yields  the  optimal  frequency- 
dependent  scales  for  the  given  suboptimal  con¬ 
stant  scale.  The  robustness  bound  after  each 
of  the  two  steps  is  shown  in  Figure  1.  Note 
that  the  peak  of  the  bound  was  minimized  in 
the  first  step  but  at  other  frequencies  the  bound 
was  minimized  in  the  second  step.  The  optimal 
frequency-dependent  scales  dAi(u;)  and  dA2i<^) 
obtained  after  the  second  step  are  shown  in  Fig¬ 
ure  2.  The  constant  scale  Do  is 
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1.0000  -0.0027 

-0.0027  0.2094 


Figure  1:  Example  1,  DD  Robustness  Bounds 


and 


A  = 


A^nc  0 

0  ^perf 


Da{u}) 


dAi(w)  0 
0  dA2{<^)I  ’ 


The  robustness  bound  obtained  using  the  DD 
procedure  can  be  compared  with  the  limiting 
cases  obtained  by  treating  all  uncertainty  as  TV 
and  using  optimal  constant  scales  and  by  treat¬ 
ing  all  uncertainty  as  LTI  and  using  optimal  fre¬ 
quency  dependant  scales.  These  limiting  cases 
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Figure  2:  Example  1,  Frequency  Dependent 
Scales 

are  shown  in  Figure  3.  The  bound  obtained  us¬ 
ing  the  DD  analysis  technique  is  between  the 
other  two  bounds  at  each  frequency  as  it  should 
be.  In  this  example,  the  DD  bound  is  essentially 
equal  to  the  bound  obtained  using  all  frequency 
dependant  scales.  It  can  be  seen  that  assuming 
all  TV  uncertainty  would  be  conservative. 


Figure  3:  Example  1,  Robustness  Bounds  for 
Different  Classes  of  Uncertainty 


4.2  Example  2 

In  this  example  the  closed  loop  system  obtained 
after  the  fourth  DKD  iteration  in  [9]  is  analyzed 


and  compared  to  the  results  presented  there.  The 
robustness  bound  obtained  for  this  system  after 
the  first  and  second  steps  described  in  Section  3 
are  shown  in  Figure  4.  As  in  the  first  example, 
it  can  be  seen  that  the  peak  of  the  bound  is  not 
lowered  in  the  second  step.  The  DD  analysis 
results  are  compared  with  the  limiting  cases  as 
discussed  in  the  previous  example  and  are  shown 
in  Figure  5.  The  analysis  of  this  system,  which 
has  both  TV  and  LTI  uncertainty,  demonstrates 
that  treating  all  uncertainty  as  TV  can  be  con¬ 
servative,  while  treating  it  all  as  LTI  can  give 
misleading  robustness  bounds.  Thus,  in  order 
to  get  accurate  analysis  information,  the  type  of 
each  uncertainty  must  be  properly  accounted  for 
by  the  analysis  procedure. 


Figure  4:  Example  2,  DD  Robustness  Bounds 

In  [9],  the  DKD  procedure  outlined  in  Sec¬ 
tion  2  was  used  for  analysis.  The  robustness 
bound  using  this  technique  is  shown  in  Figure  6 
along  with  the  robustness  bound  obtained  using 
the  DD  analysis  procedure.  The  peak  of  the  DD 
bound  is  less  than  the  peak  of  the  bound  in  [9] 
because  the  constant  and  frequency-dependent 
scales  that  minimize  the  peak  of  the  robustness 
bound  were  found  simultaneously,  whereas  in  the 
DKD  procedure,  the  constant  scales  are  found 
independently  of  the  frequency-dependent  scales. 
When  the  DKD  procedure  is  used  in  an  itera¬ 
tive  manner  as  in  [9],  the  constant  scales  found  in 


7 


1.5 


c 

D> 

ns 


0.5h 


Constant 
I  —  Mixed 

Frequency  Dependent 


10 


-2 


10  10  10 
Frequency  (rad/sec) 


10 


Figure  5:  Example  2,  Robustness  Bounds  for 
Different  Classes  of  Uncertainty 


later  iterations  are  improved  due  to  the  frequency- 
dependent  scales  which  have  been  absorbed  into 
the  system.  Evidence  of  this  can  be  seen  in  Fig¬ 
ure  7,  where  the  constant  scales  are  found  with¬ 
out  the  aid  of  frequency-dependent  scales  from 
a  previous  iteration.  However,  even  when  the 
DKD  procedure  is  used  iteratively  until  it  con¬ 
verges  as  in  [9]  the  results  can  be  conservative  as 
shown  here. 


Figure  6:  Example  2,  Robustness  Bounds  using 
DKD  and  DD 

We  also  compare  the  DD  robustness  bound 
with  the  bound  obtained  using  the  ZDs  proce¬ 
dure  outlined  in  Section  2.  The  ZDs  procedure  is 


the  only  feasible  analysis  alternative  to  the  DD 
procedure  for  general  systems  containing  both 
TV  and  LTI  uncertainty,  that  is,  for  systems 
for  which  the  controller  was  not  synthesized  us¬ 
ing  the  technique  of  [2]  or  [9].  The  bounds  are 
shown  in  Figure  7.  Again,  the  peak  of  the  ro¬ 
bustness  bound  obtained  by  the  DD  procedure 
is  less  than  that  obtained  by  an  alternative  pro¬ 
cedure.  Here  the  difference  is  even  more  dra¬ 
matic  than  shown  in  Figure  5.  This  is  because 
as  with  the  DKD  method,  the  ZDs  procedure 
finds  the  constant  scales  independent  of  the  fre¬ 
quency  dependent  ones,  but  unlike  the  DKD 
procedure,  the  ZDs  method  doesn’t  benefit  from 
frequency-dependent  scales  obtained  in  a  previ¬ 
ous  iteration.  It  can  be  seen  that  the  DD  anal¬ 
ysis  procedure  can  remove  a  significant  amount 
of  conservatism  when  analyzing  a  general  system 
with  both  TV  and  LTI  uncertainty. 


Figure  7:  Example  2,  Robustness  Bounds  using 
ZDs  and  DD 


5  Summary 

A  new  technique  for  computing  the  robustness 
bound  for  systems  with  both  TV  and  LTI  un¬ 
certainties  was  presented.  The  method  provides 
less  conservative  results  than  current  techniques. 
Furthermore,  it  can  be  used  on  any  system  that 
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has  TV  and  LTI  uncertainty  in  linear  fractional 
form. 

While  the  numerical  results  presented  consid¬ 
ered  real  scales,  complex  scales  could  have  been 
used  as  well.  However,  in  [6]  and  [10]  it  was 
shown  that  allowing  the  constant  scales  to  be 
complex  doesn’t  reduce  the  scaled  singular  value 
if  the  system  being  analyzed  has  only  real  state 
space  entries.  Furthermore,  in  [11],  Packard  shows 
that  restricting  scales  associated  with  full  com¬ 
plex  LTI  blocks  to  be  positive  real  can  be  done 
without  loss  of  generality.  That  leaves  only  the 
case  of  repeated  LTI  uncertainty  with  hope  of 
benefitting  from  the  use  of  complex  scales.  There 
may  be  cases  where  this  improvement  is  signifi¬ 
cant  but  it  hasn’t  been  a  factor  in  the  problems 
we  have  encountered. 

Since  the  first  step  of  the  analysis  is  more 
computationally  intensive  than  the  second  it  is 
often  advantageous  to  obtain  the  suboptimal  con¬ 
stant  scales  using  a  few  strategically  selected  fre¬ 
quency  points.  Then,  step  two  can  be  performed 
using  a  finer  frequency  grid  to  give  a  robustness 
bound.  It  should  be  noted  that  using  this  ap¬ 
proach  it  is  possible  for  the  peak  of  the  bound 
obtained  in  the  second  step  to  be  larger  than  the 
peak  obtained  in  the  first  step,  thus  indicating  a 
frequency  point  that  should  be  included  in  the 
first  step. 

Finally,  while  a  set  of  LMIs  is  mathematically 
equivalent  to  a  single  LMI  containing  the  set  of 
LMIs  on  its  diagonal,  we  have  found  that  the 
first  case  requires  significantly  less  memory  and 
often  is  much  faster  to  solve. 
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A  Appendix 


The  state  space  matrices  A,  B,  C,  and  D  of  the 
weighted  closed-loop  system,  G,  used  in  exam¬ 
ple  one  are  given  below,  where  the  closed  loop 
matrix  F„(G,  A)  =  A  -t-  BA{I  —  DA)~''-C. 
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